Far-Field Electron Spectroscopy of Nanoparticles 



O 

o 

(N 
< 

in 

(N 



c3 



13 

c 

o 

(N 
> 

in 

(N 

m 
O 

o 

-i— > 

c3 



i 

C 

o 
o 



X 



M. A. Itskovsky 1 , H. Cohen 2 and T. Maniv 1 
1 Schuhch Faculty of Chemistry, Technion-IIT, 32000 Haifa, ISRAEL 
2 Weizmann Institute of Science, Chemical Research Support, Rehovot 76100, ISRAEL 

(Dated: February 5, 2008) 

A new excitation mechanism of nanoparticles by relativistic, highly focused electron beams, is 
predicted under nanoparticle recoil. The corresponding electron energy loss spectra, calculated for 
metallic (silver and gold) and insulating (Si02) nanoplatelets, reveal dramatic enhancements of 
radiative electromagnetic modes within the light cone, allowed by the breakdown of momentum 
conservation in the inelastic scattering event. These modes can be accessed with e-beams in the 
vacuum far-field zone, offering interesting capabilities for signal transmission across nano to meso 
length scales. 

PACS number(s): 79.20.Uv, 78.67.Bf, 73.20.Mf, 41.60.-m 



A powerful approach to the study of nanoparticles is 
provided by very fast (relativistic) electron beams (e- 
beams), with typical lateral resolution on an atomic scale, 
available in scanningtransmission electron microscopes 

(stem) ma a iiaB- 

As discussed previously [H, [3,13] , 
when the e-beam is restricted to the vacuum near a 
selected nanoparticle, its electromagnetic (EM) interac- 
tion with surface plasmons or surface plasmon-polaritons 
(SPPs) is reminiscent of the near-field interaction of 
subwavelength optical probes. Several works have re- 
cently studied realizations of Cherenkov radiation exci- 
tation within various dielectric media by e-beams moving 
in near-field vacuum zones @, 0, [l^, 11 1 . 



In the present paper we show that by considering re- 
coil effects of the nanoparticle during the scattering event 
one introduces far-field coupling between the electron and 
the nanoparticle, which dramatically enhances the radia- 
tive channels in the loss spectrum. To illustrate our main 
points we consider here a simple model where the e-beam 
is propagated in the vacuum along a wide face of a rect- 
angular nanoplatelet (oriented, e.g., in the x — y plane), 
and a surface or guided wave induced by the electron 
is propagated with a wave number k x along the beam 
axis. The spatially sensitive nature of the correspond- 
ing electron energy loss process arises from the exponen- 
tial dependence, e~ 2K b , of the EM interaction between 
the e-beam and the platelet on the impact parameter 

b. The extinction coefficient, K* 



\jk 2 — (lu/c) 2 , with 



k 2 = k 2 + k 2 , determines the tail of the evanescent field 
in the vacuum for values of k outside the light-cone, i.e. 
for k > lu/c. Inside the light-cone, i.e. for k < lu/c , 
K* is purely imaginary and the corresponding interac- 
tion becomes spatially oscillating, allowing the electron 
to exchange photons with the particle far away into the 
vacuum. This striking possibility has been overlooked 
in the recent literature of STEM-electron energy loss 
spectroscopy (EELS), since the excitation by an electron 
moving in the vacuum with a classical velocity v, has 
been restricted to a constant longitudinal wavenumbcr 
k x = lu/v > cj/c, implying EM coupling to the nanopar- 



ticle which is restricted to the evanescent tail near the 
surface. 

Following Ref . [4] , the focused e-beam is described here 
as a one-dimensional wave, propagating along the x-axis, 
while in the transverse (y — z) directions it is described 
by a wave-packet localized within a smoothly converg- 
ing cross section along the beam axis, whose shape is 
assumed to be squared for the sake of simplicity. The 
nanoplatelet half-sides a*, b*, c* along the x, y, z axes re- 
spectively, are assumed to satisfy b* ^> a* ^> c* (see 
below for more details). Within the framework of this 
quantum-mechanical approach we use an effective beam- 
nanoparticle interaction Hamiltonian based on the re- 
tarded four-vector image potential of the e-beam and 
Born-Oppenheimer-like separation of coordinates (i.e. 
the 'slow' transverse coordinates from the 'fast' longi- 
tudinal one), justified by the small e-beam converging 
angle. The corresponding rate of change of inelastic scat- 
tering probability of the incident e-beam at an impact 
parameter b and for energy loss AE — Hu), calculated 
to first order in the effective interaction Hamiltonian, 
-e$ — e— A T , can be written in the form: 

mc 
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EM Vi Z > 0) 



where $ and A x are the scalar and x-component of 
the four-vector potential respectively, dominating the 
interaction with the e-beam, p x — Hq x , v — hq x /m 
its initial longitudinal velocity and m is the relativis- 

In Eq.© 



tic electron mass: m — mo/y 1 — (v/c) 
Aq x w (lu/v) + h I ql r \ — (<Zt r )" /2mv is the longitu- 
dinal momentum transfer of the e-beam, ~q\ r = [q i y ,q\ / 



2 



and ~q{ r its transverse momenta, initial and final respec- 
tively, and e lq:cX x~q tT (y, z\ x ) is an unperturbed e-beam 
eigenfunction with a longitudinal wave number q x and 
asymptotic transverse wavevector ~~qtr- Note that the 
width /3 _1 of the Gaussian distribution function, intro- 
duced in Eq. |T]) to account for the high transverse-energy 
cutoff caused to the e-beam by the objective aperture, 
determines a region around the beam focal plane whose 
length L may be used as a normalization factor for the 
e-beam wave functions in our model. 

The condition a* 3> c* ensures that the interaction 
potential between the platelet and an external electron 
at a projected distance |x| from the platelet center de- 
cays to zero at least as quickly as l/x 2 for |x| > a* 
(see, e.g., Ref. 13]). Under these circumstances the 
limits of integrations over x in Eq.(pQ) may be set at 
—a* and a*, rather than at — L and L. The poten- 
tial correlation function for t > may be expressed in 
terms of the relevant components of the 4-tensor pho- 
ton propagator D v ^(r', r ;t) , v,(i — 0,1,2,3 (<-> 

ct, x, y, z), by: -i (h p £ m (x', y', z'; t) H p j? u (x, y, z; 0) 



in the vacuum outside the rectangular platelet 0: 



nf 



£>o, ("?', r>; t) + %- c D 0A r- 1) + 2£D h o 
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(mcfi Di.i ( r ' 1 r ^ 0r ^ e sa ^ e °^ siuiphcity we ne- 
glect the finite size effects parallel to the wide (x — y) 
face of the platelet in the calculated photon propagator. 
Consequently the energy loss rate in Eq.|T|) reduces to: 



R oc dk 



dk y Im 



r {k x , k y ,u) 
K* 



e'^) 2 ' 2m \l (yl, T tr ; ky,K*; Aq x - k, 



1 tr< 1 ir 

where: 



(2) 



il /l dXel(AqX ~ kx)XJ ( ^*ri ^ K *' X ) ' ^ 



and Jyqtn qU k y> K *'i x ) 

f dze K * z f dye~ lk yy X *-?f {v,z;x)x-?* (y,z;x). The 

y tr 

effective surface dielectric response function, appearing 
in Eq.Q, is given by: 

r(k x ,k y ,co) « r ,o + (hql/mc) r ,i + 

(hqt/mc) n, + (ft 2 qlq x / (mc) 2 J n,i 

where r„ iM are the components of the EM reflection 4- 
tensor, which determine the dressed photon propagator 



5 v „e- K '\ z -*'\-r u „e K *(*+ z ') 



with 77,, = 1 , or —1 for !/ = , or v = 1, 2, 3 respectively. 
In the long wavelengths limit discussed in Ref. 12j we find 
that: 



Im 



r (lc,w) /K* 



Im 



{ \{K*/k 2 ) U + {{v/cf - (co/ck) 2 ) f /K* 



-2K* 



(4) 

where f e = (e 2 K* 2 - Q 2 ) / D+ D~ , f = 
(K* 2 - Q 2 ) /D+D- , D+ = eK* + Q tanh (Qc*) , D~ = 
eK* + Qcoth(Qc*) , L>+ = AT* + Qtanh (g_c *) , D~ = 
K* + Qcoth (Qc*) , Q = ^fc 2 - (J/c) 2 eH, and e (w) 
is the local bulk dielectric function of the platelet. In 
the limit of a semi-infinite medium the resulting expres- 
sion reduces to the surface dielectric response function 
obtained in Ref. [3] by using Maxwell's equations with 
macroscopic boundary conditions. 

The standard classical approximation for the loss 
function[3| is obtained from Eq.© by making the fol- 
lowing assumptions: (1) the e-beam transverse momen- 
tum distribution function J ( ~~q{ r , ~~ql r ', k y , K*; x\ is a 

constant, that is equivalent to a 5-function in the corre- 
sponding real-space transverse coordinates, (2) the con- 
tribution of the transverse energy to Aq x (see below 
Eq.([T])) can be neglected, and (3) the effective parti- 
cle size, a*, appearing as an integration limit along the 
beam axis, is infinite. A simple estimate shows, how- 
ever, that assumption (2) is usually invalid since the 
contribution of the transverse energy to Aq x can be as 
large as (oj/v). Assumption (3), in conjunction with (1), 
yields the conservation of longitudinal momentum, i.e. 
Aq x — k x = 0, which together with assumption (2) im- 
poses the fixed condition k x = (lo/v). In the present 
paper we remove only the third assumption by allow- 
ing a* to be a finite length, which reflects an effective 
range of the beam-particle interaction along the beam 
axis. Consequently the longitudinal momentum distri- 
bution around Aq x — k x =0, defined by the integral in 
Eq.([3]), is smeared and many wavenumbers k x inside the 
light-cone start contributing to the loss rate, Eq.©. 

The condition for the smearing to be significant is 
7r / 'a* > (lo/v), so that typically for frequencies to in 
the visible range a* should be smaller than 200 nm. 
Nanoplatelets of that lengths should dramatically en- 
hance radiative excitations by the e-beam, previously 
overlooked in the literature. It should be stressed that, 
for the sake of simplicity, this is done by invoking the 
dielectric response function of an infinitely wide platelet. 
A consistent treatment of the breakdown of translation 
invariance is expected, however, to further enhance all 
radiative channels. 
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FIG. 1: EEL spectra (solid lines) of 100 nm long Ag and Au 
platelets calculated for a 100 keV external e-beam at various 
impact parameters between 10 and 40 nm. The experimen- 
tal optical dielectric functions, e(u>), for silver and gold |15j ] 
have been exploited. Dashed lines represent spectra calcu- 
lated by the classical theory. Inset: SPP dispersion curves, 
uj (Re k) , uj (Im k) in the complex fc-plane for silver. The in- 
dicated values of k and w are normalized by k n = u n /c, and 
uj n = 10eV, respectively. 



As a first example we calculate the EEL function of a 
100 nm long silver and gold platelets for an external 100 
keV e-beam at various impact parameters (see Fig.l). 
To analyze the various SPP resonances one may consider 
the zeros of the denominator of the extraordinary wave 
amplitude f e in Eq.(0J in the complex fc-plane. With 
the experimental optical dielectric function, s(lo), for sil- 
ver 



15j the resulting dispersion relation (inset, Fig.(l)) 



exhibits a rather flat branch of cj(Rcfc) inside the light- 
cone, which can be attributed to radiative SPP, seen as 
a mirror image of the usual non-radiative SP dispersion 
curve with respect to the light-line. The sector of w(Re k) 
connecting the two branches across the light line exhibits 
a negative slope, where Imfc(w) oc Imefw) has a sharp 
peak. The sharp dip in the EEL spectrum just above the 
classical SP frequency (at 3.8 eV) reflects these closely 
related features. These peculiar features are missing in 
the loss spectrum of the gold platelet, Fig.l. 

At slightly higher frequencies the EEL signals exhibit 
a pronounced rise due to the enhanced SPP density of 
states associated with the flat radiative SPP branch. The 
corresponding EEL peak intensity exhibits a remarkably 
weak attenuation, in marked contrast to the exponential 
attenuation of the main SP peak, calculated in the classi- 
cal limit for increasing impact parameter. The radiative 
nature of the beam-particle coupling shown in Fig.(l) is 
even more pronounced in the low energy region below 



FIG. 2: EEL spectra (solid lines) of a 100 keV e-beam prop- 
agating parallel to the x-axis of a rectangular Si02 platelet 
at distances b — 10, 20, 30 nm above its wide [x — y) face. 
The platelet half sides along the x, and z axes are: a* — 50 
nm, and c* = 5 nm respectively. The corresponding spec- 
tra (dashed lines) obtained from the classical theory are also 
shown for comparison. Note the close similarity of the classi- 
cal spectrum for b — 10 nm to that obtained in Ref. [lG] for the 
same e-beam velocity at nearly the same impact parameter 
parallel to a sharp wedge (see inset). Inset: A schematic illus- 
tration of the Si02 platelet showing the e-beam (wide face) 
configuration studied in the present paper and the (wedge) 
configuration investigated in Ref. [3]. 



the main SP peak, where the classically calculated signal 
drops to zero. Here our calculated EEL function exhibits 
a pronounced broad band with linearly increasing inten- 
sity for increasing frequency and almost no attenuation 
with increasing impact parameter. These features are 
due to the fact that the loss signal well below the main SP 
frequency is dominated by the contribution from the or- 
dinary wave amplitude / Q , appearing in Eq.(|3]), which is 
singularly enhanced near the light line (where K* — > ) , 
and thus reflecting the nearly pure (transverse) photonic 
nature of the excitations by the e-beam in this 'classically 
forbidden' region. 

A similar and even more dramatic situation arises in 
the forbidden energy gap of the electron-hole pair excita- 
tions in semiconductors and insulators, as shown, e.g. in 
Fig. (2) for a 100 nm long Si0 2 platelet. The EEL spec- 
tra of an external 100 keV e-beam, propagating parallel 
to the x — y face of the platelet at various impact pa- 
rameters b, reveal a pronounced broad peak within the 
forbidden gap region, which does not decay with increas- 
ing b values. Similarly to the situation with the silver 
platelet well below the main SP peak, the strong radia- 
tive nature of this feature arises from the ordinary wave 
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amplitude f a , corresponding to the excitation of purely 
transverse EM waves, polarized within the x — y plane, 
which totally dominates the loss signal in the gap region. 

The spectrum shown in Fig. (2) for an impact param- 
eter b = 10 nm may be compared to the result reported 
in Ref.[lj| for an electron moving parallel to a 90° SiC>2 
wedge at a distance of 8.5 nm (see Inset to Fig. (2)). The 
pronounced radiative broad band within the gap region, 
obtained in our calculation, dramatically contrasts the 
vanishing loss signal shown there in Fig. (4) for an elec- 
tron with the same velocity ( v = 0.54c) and nearly 
the same impact parameter. The lack of far-field cou- 
pling in the latter theoretical approach restricted the 
fast external e-beam to excite EM waves propagating 
within the dielectric medium onlyQ, similar to ordinary 
waveguide modes which can develop within a thin Si02 
slab in the forbidden gap region where Re e (u>) ~ 2 , 
and Ime(w) — * 0. For ideal planar geometry (as as- 
sumed in our calculation of the dielectric response func- 
tion r ^ k ,Luj), the corresponding waveguide modes ap- 
pear as extremely narrow resonances which can not be ex- 
cited by an e-beam with Aq x values outside the light cone 
due to the vanishingly small dielectric damping, Ime (a;). 

Such radiation excitations become possible for the non- 
planar geometries studied in Refs.0,[l3| even under the 
recoilless scattering approximation exploited there (but 
only above a threshold beam energy considerably higher 
than 100 keV) due to the translational symmetry-broken 
dielectric media considered in their calculations. Yet, the 
corresponding Cherenkov channels remain fundamentally 
different from the ones we propose: The opening of scat- 
tering channels with wave numbers inside the light cone 
allows coupling of the e-beam to the continuum of EM 
modes which are extended into the vacuum, leading to 
dramatic enhancement of the loss signal at e-beam en- 
ergies around 100 keV. The relative strength of the ra- 
diative effect calculated here may be further appreciated 
by noting the close proximity of the e-beam to the sharp 
SiC>2 wedge in [16j .which strongly enhances the (evanes- 
cent) near-field coupling with the e-beam. 

In summary, applying a quantum-mechanical approach 
to the scattering problem, we have shown that Cherenkov 
radiation of highly focused relativistic e-beams in STEM, 
discussed recently in the literature [1, 0, [l(|, have a much 
broader scope than originally presented. The dramatic 
enhancement of radiative channels, imposed by the finite 
nanoparticle size, arises from the breakdown of momen- 
tum conservation in the inelastic scattering event along 
the e-beam axis. Further enhancement should be realized 
due to the extreme lateral confinement of the e-beam and 
its associated transverse momentum uncertainty, which 
were neglected here. The radiation predicted to be emit- 
ted from both conducting and insulating nanoplatelets 



can be generated at impact parameters much larger than 
the evanescent tail of the excited surface EM modes. 
Large deviations from the classical, nonradiative EEL 
signal are found to persist also at small impact param- 
eters, which can be readily tested experimentally. The 
proposed spectroscopic experiments are closely related 
to recent developments in the field of SPP-based far-field 
optics. 17|. Far field EELS can therefore become a use- 
ful complementary tool for an already recognized new 
emerging field. 
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